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SECTION  1 


INTRODUCTION 


This  is  the  third  in  a  series  of  reports^’  *  documenting  a  useful,  approxi¬ 
mate  solution  of  the  Boltzmann  Equation  for  the  energy  spectrum  of  swarm  electrons 
(energy  less  than  a  few  tens  of  eV)  in  a  weakly  ionized  gas.  We  review  the  theory 
and  then  discuss  time-dependent  problems,  extending  the  steady-state  discussions  of 
References  1  and  2. 

The  electrons  are  t<iken  to  suffer  realistic  momentum-transfer  collisions  (due 
to  both  elastic  and  inelastic  scatterings),  «ind  inelastic  energy-transfer  collisions  with 
the  background  host  gas  molecules.  Requisite  cross  sections  are  taken  from  data. 
Electron-electron  and  electron-ion  collisions  are  ignored  compared  to  electron-neutral 
collisions.  The  gas  molecules  are  in  their  own  thermal  equilibrium;  rotational  levels  of 
diatomic  molecules,  for  example,  are  thermally  populated  with  a  Maxwell-Boltzmann 
distribution  at  ambient  temperature,  and  super-elastic  collisions  of  electrons  with 
excited  states  are  included.  We  do  not  include  thermally  populated  vibrational  states 
because  of  unknown  cross  sections,  although  in  principle  they  could  easily  be  included. 

An  electric  field  E  is  present  which  causes  the  electrons  to  drift  and  to  absorb 
energy  from  the  field.  It  makes  the  distribution  function  non-isotropic.  If  the  electric 
field  is  constant  in  time,  the  energy  gained  from  it  is  lost  to  energy-transfer  collisions 
at  the  same  average  rate,  so  the  swarm  will  awrhieve  a  steady-state,  which,  however,  is 
not  a  Maxwellian.  The  gas  is  considered  an  infinite  energy  sink. 

Reference  1  presented  the  basic  formulation,  and  introduced  the  usual  mean 
Loss  Function  L,  and  the  Straggling  Function  >1.  It  was  shown  that  the  energy  spec¬ 
trum,  and  most  swarm  transport  coefficients,  are  determined  to  a  good  approximation 
by  the  momentum  transfer  cross  section  and  two  moments,  L  and  Af,  of  the  inelastic 
energy  transfer  cross  sections.  The  spectrum  is  reduced  to  quadratures  over  these 
three  molecular  parameters.  Reference  1  discussed  the  physical  importance  of  strag¬ 
gling,  obtained  the  differential  equation  for  the  time  evolution  of  the  energy  spectrum 
in  time  varying  fields,  presented  the  expression  for  the  steady  state  spectrum  in  a 
constant  field,  zmd  applied  it  to  N2,  comparing  the  spectrum  itself  with  more  detailed 
calculations,  and  comparing  computed  swarm  parameters  with  data.  Agreement  with 
most  swarm  parameters  is  striking  over  some  five  orders  of  magnitude  in  the  experi¬ 
mental  parameter  E/Nq.  (No  is  the  molecular  number  density.)  An  exception  is  the 


avalanche  rate,  which  is  rather  larger  than  data  at  relatively  small  E/Nq.  (Accurate 
agreement  with  the  avalanche  rate  would  not  be  expected.) 

The  form  for  the  steady  state  spectrum  is  easily  understandable  in  terms  of 
ordinary  diffusion-convection  theory,  being  essentially  a  Fokher-Planck  expansion  in 
energy  space,  and  requires  small  fractional  mean  energy  loss  per  collision.  Fractional 
energy  loss  is  sufficiently  small  in  Nj  and  especially  O2  for  the  approximation  to  be 
valid  over  a  useful  wide  span  of  E/Nq.  Reference  1  also  explained  the  physical  reason 
for  the  inadequacy  of  the  Continuous  Slowing  Down  Approximation  (CSDA)  [in  which, 
for  example,  the  avalanche  rate  is  orders  of  magnitude  too  small  at  relatively  small 
E/Nol. 

After  a  brief  historical  sketch.  Reference  2  showed  that  the  approximate 
steady  state  spectrum  becomes  exact  in  the  limit  that  molecular  quantum  transition 
energies  are  small  compared  with  the  swarm  energy  and  the  gas  temperature.  The 
proof  relies  on  detailed  balance.  This  further  shows  that,  in  addition  to  the  physical 
reasons  for  its  inadequacy,  the  CSDA,  which  is  commonly  employed  in  analytic  treat¬ 
ments,  suffers  the  formal  deficiency  that  it  violates  detailed  balance.  The  steady  state 
spectrum  was  then  used  to  compute  Oj  swarm  parameters,  and  compared  with  data. 
Again  the  agreement  is  striking  over  some  five  orders  of  magnitude  in  E/Nq,  with  the 
avalanche  rate  again  being  an  exception.  In  addition,  an  expression  was  presented  for 
the  energy  spectrum  in  a  sinusoidally  time  varying  electric  field.  Finally,  some  more 
fundamental  issues  were  raised  concerning  basic  statistical  mechanical  approaches  to 
the  problem  of  determining  the  electron  energy  spectrum  under  the  realistic  random 
walk  behavior  in  energy  space. 

Recommendations  for  some  needed  cross  sections  were  made  in  References  3 

and  4. 

In  the  present  report,  we  first  review  the  basic  theory  in  Section  2.  Then  in 
Section  3  we  present  a  numerical  finite  difference  algorithm  for  solving  the  differen¬ 
tial  equation  previously  derived  as  an  approximation  to  the  time-dependent  integro- 
differential  Boltzmann  Equation.  The  remaining  sections  present  veu-ious  numerical 
solutions  for  time-evolving  spectra. 
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SECTION  2 


BASIC  EQUATIONS  AND  REVIEW  OF  THEORY 


In  this  section  we  briefly  review  the  theory  for  the  approximate  Boltzmann 
Equation  and  its  steady  state  solution.  Details  are  in  References  1  and  2. 


#  2.1  THE  BOLTZMANN  EQUATION 


Dropping  spatial  gradient  terms,  the  Boltzmann  Equation  for  the  electron 
distribution  function  f(v,t)  in  an  electric  field  E  and  in  a  background  gas  is 


at 


—  •  Vyf  =  C  +  Sources  —  Sinks 
m 


(!)• 


Here  e  is  the  magnitude  of  the  electron  charge.  Actual  sources  amd  sinks  are  problem 

•  specific.  They  do  not  play  a  central  role  here,  and  will  be  dropped  until  a  specific 
problem  is  considered.  We  concentrate  on  the  effects  of  E  and  the  collision  integral  C 
in  determining  f. 

Many  bulk  plasma  parameters  depend  only  on  the  speed  distribution  func- 

^  tion 

F(v)  =  f  dn,f(vl  (2) 

#  or,  equivalently,  on  the  energy  spectrum  g(w). 


g(w)  =  ^F(v) 
m 


(3) 


#  where  w  =  mv*/2  is  the  electron  energy.  If  n 

normalizations  are^ 


is  the  electron  number  density,  the 


n  =  f(^d\  =  j  F(v)v*dv  =  j  gdw 


(4) 


3 


By  suitable  approximation  to  the  collision  integral,  we  shall  obtain  a  differential  ^ 

equation  for  F  or  g.  F  will  also  be  considered  a  function  of  w  normalized  to 
/  F(w).y/wdw  =  n. 


2.2  LEGENDRE  EXPANSION  FOR  f  • 

For  reduced  electric  fields  E/Nq  not  too  large,  Vdrift  <  Vthermah  and  the 
departure  from  isotropy  will  be  small.  Then  f  is  usefully  developed  in  the  usual 
Legendre  Polynomial  series  in  the  angle  0  between  v  and  E  [References  6,7].  ^ 


f  = 

t 

—  fo  +  fi  cos  0  -h  small  terms  (5) 

When  electron-neutral  collisions  dominate,  C  is  linear  in  f.  C  is  also  expanded, 
C  =  Co  +  Cl  cosd  4-  •  •  •,  and  when  the  series  is  inserted  in  Eq.  (5),  the  Po  and  Pi 
parts  become®'^ 


afo 

dt 


^  (2^-i  +  fill 

3m  V  V  av  J 


—  Co 


afi  eE  ^0 
dt  m  dv 


fo  is  essentially  the  energy  spectrum,  for 

F(v)  =  j  {io  +  ficos0)dnv  =  47rfo 
while  fi  controls  bulk  electron  motion  with  drift  velocity 


Vd 


-  I  fvd®v 
n  J 


-  j {io  +  fi  cos  0)vd^v 


(6a) 

(6b)  • 


(7) 


(8) 
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where  E  is  a  unit  vector  in  the  direction  of  E. 


2.3  LEGENDRE  EXPANSION  OF  COLLISION  INTEGRAL 


Thus  Co  is  the  collisional  contribution  to  the  time  rate  of  change  of  the 
energy  spectrum.  It  is  discussed  in  Section  2.6  below. 

Cl  is  the  contribution  to  the  rate  of  change  of  that  part  of  the  spectrum  that 
controls  the  swarm  momentum.  In  general,  for  a  weakly  ionized  gas  in  which  electron- 
neutral  collisions  dominate®’^,  and  when  elastic  collisions  dominate  the  momentum 
transfer, 


Cl  =  (9) 

where  Um  =  Nov^Tm  is  the  mono-energetic  momentum  transfer  collision  frequency, 
is  the  momentum  transfer  cross-section,  and  No  is  the  background  gzis  number  density. 


♦  2.4  SLOWLY  VARYING  FIELDS 


Equations  (6)  with  (9)  and  a  model  for  Co  could  be  solved  together  forward 
in  time  to  obtain  the  evolution  of  the  distribution  function,  and  so  all  transport  coef¬ 
ficients.  However  a  single  equation  can  be  obtained  when  fields  vary  slowly  compared 
with  Urn-  Inserting  Eq.  (9)  into  (6b), 

^1  eEdfo  _ 

o.  o  ^mfl 

at  m  av 

Thus,  when  d/dt  C  Um,  the  first  order  distribution  function  is  expressible  in  terms  of 
the  instantaneous  energy  spectrum  without  any  time  lag: 

f  - 

‘  rni/m  dv 

Since  is  generally  greater  than  some  10"sec“'  in  sea  level  air,  condition  Eq.  (11) 
holds  in  many  circumstances. 
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The  single  equation  for  fo  is  now  obtained  by  combining  Eq.  (11)  with  (6a). 
After  a  little  algebra,  it  becomes 


dt  3m*v*  dv  dv )  ° 

This  may  be  written  as  well  as  an  equation  for  g(w), 


(12) 


^  A  /f!^ 

dt  dw  (mi/rn 


(13) 


where  Cg  is  the  collision  integral  for  g.  In  this  equation  we  have  split  the  electric 
field  contribution  to  the  energy  flux  (curly  brackets  in  Eq.  (13))  into  two  parts,  one 
proportional  to  g,  and  one  to  5(wg)/5w.  Equation  (13)  exhibits  the  field’s  role  in 
determining  g. 


2.5  ENERGY  DIFFUSION  IN  THE  ELECTRIC  FIELD 


The  two  parts  correspond  to  a  mean  flux  (^g),  and  a  diffusive  flux 
(~  5(wg)/5w).  As  is  well  known,  elastic  scattering  in  the  presence  of  an  electric 
field  produces  both  a  mean  gain  in  energy  (rate  ~  e^E^/rni/m),  and  a  diffusion  in 
energy  (diffusion  coefficient  ~  (2/3)we*E*/mi/m)-  Diffusion  is  caused  by  scatterings 
parallel  or  anti-parallel  to  E,  spreading  out  the  electrons’  positions  along  the  field, 
and  therefore  their  potential  energy. 

That  the  diffusion  term  is  important  can  be  seen  by  comparing  its  effect  in 
Eq.  (13)  to  that  of  mean  energy  gain  on  a  Maxwellian  distribution  of  temperature  T: 

Diffusive  Gain  Rate  _  —  |^(wg)  _  2  w 

Mean  Gain  Rate  g  ^  3  T 

This  is  of  o  -li-x  unity,  and  is  therefore  important.  It  changes  sign  over  the  average 
energy  3T/2.  T-  'sion  allows  pau^ticles  of  energy  w  >  3T/2  to  gain  energy  faster 
than  the  i  ean,  ;..id  particles  with  w  <  3T/2  to  gain  energy  slower  than  mean,  thus 
spreading  out  the  spectrum. 
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1 

2.6  THE  INELASTIC  COLLISION  INTEGRAL 

The  general  form  for  the  collision  integral  is 

Cg  =  y*  g(w')P(w',  w)dw' —  g(w)  j  P(w,  w')dw' 

(15) 

• 

where  P(w,w')dw'  is  the  probability  per  second  that  an  electron  with  energy 
and  ends  up  with  energy  (w',w'  +  dw').  Inelastic  collisional  excitations 
molecular  states  k  of  energy  w^  and  cross  section  tyk(w)  contribute  to  ? 

w  collides 
to  discrete 

an  amount 

P(w,w')dw'  =  NoV^ffk5(w  —  Wk  —  w')dw' 

k 

(16) 

• 

assuming  all  host  molecules  in  the  ground  state.  (This  assumption  will  be  relaxed  later. 

Rotational  levels  will  be  thermally  populated.)  Inserting  this  in  Eq.  (15)  produces  the 
usual  expression, 

€ 

Cg  =  No  Xll(^kVg)w+,^  -  (^kVg),) 
k 

(17) 

m 

Since  the  cross  sections  are  not  exact  step  functions,  the  sum  Cg  itself  is,  of  course,  a 
continuous  function  of  w.  Developing  the  summand  in  a  Taylor  Series  about  w,  we 
exhibit  the  first  two  terms. 

(18) 

• 

where 

L  =  X)wkak(w) 
k 

(19a) 

• 

M  =  X]  ^wjak(w) 

k  ^ 

(19b) 

• 

L  is  the  Loss  Function  describing  the  mean  energy  loss  rate  to  inelastic 
A(  is  a  measure  of  the  spread  of  energy  loss  rates  about  the  mean.  It 

7 

processes. 

is  referred 

• 

to  as  the  Straggling  Function,  and  approximates  the  effects  of  the  random  walk  in 
energy  space.  Energy  loss  to  heavy  particle  recoil  from  elastic  scatterings  can  easily 
be  incorporated  in  the  definitions.  L  and  M  are  constructed  from  measured  cross 
sections.  The  Continuous  Slowing  Down  Approximation  (CSDA)  is  that  only  the  L 
term  is  retained. 

Using  Eq.  (18),  the  equation  for  <9g/<9t  becomes 


^  4.  A 

di  dvi 


In  flux  conservative  form, 


fe^E^ 

2d,  ■ 

]  d 

d  , 

g - (wgi 

®  3  3w^ 

VgL  +  — (vgAaj 

CPN 

(20) 

• 

where 


dg  d4> 

—  +  —  =  0  , 

dt  dvf 


e^E^ 

mi/rn 


2  3.  . 

e-5g;;(wg) 


-  NovgL  -  No^(vgM) 


(21) 


(22) 


is  the  flux  in  energy  space. 

The  effect  of  thermal  motion  of  gas  scatterers  was  discussed  in  [Reference  l], 
where  it  was  shown  that  the  Davydov  term  that  usually  accounts  for  it  can  be  neglected 
if  thermally  populated  rotational  states  are  included. 


2.7  SPREAD  IN  ENERGY  LOSS 


The  ratio  M/L  is  more  slowly  varying  than  L  or  Af  itself,  since  irregular 
behavior  in  the  cross  sections,  such  «is  the  N2  resonance  spike  near  2.5  eV,  will  occur 
in  both  L  and  Af  and  divide  out  in  the  ratio. 

Spread  about  the  mean  energy  loss  rate  L  can  be  compared  with  L.  From 
Eq.  (  22)  it  is  of  order 


Diffusive  Loss  Rate 
Mean  Loss  Rate 


1  d 
vgL  3w 


(vgX)  ~  ^ 


(23) 
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where  w  is  the  energy  over  which  Mg  =  (M/L)Lg  changes  much,  being  of  order  the 
swarm  temperature  or  a  characteristic  energy  over  which  L  itself  changes.  This  ratio 
is  not  negligibly  small,  and  can  be  of  order  unity  (or  larger,  near  thermal  energies 
where  L  changes  sign),  showing  that  diffusive  contributions  to  energy  loss  can  be 
comparable  to  the  mean  loss,  as  was  the  case  for  energy  gain  in  the  electric  field. 
Since  in  steady  state  mean  loss  essentially  balances  mean  gain,  diffusive  loss  can  be 
as  important  as  diffusive  gain  in  determining  the  spectrum. 


2.8  IMPORTANCE  OF  ENERGY  DIFFUSION 


The  steady  state  spectrum  is  that  for  which  =  0,  and  results  from  a 
balance  between  total  energy  gain  from  E  and  total  energy  loss  to  collisions. 

It  is  instructive  to  temporarily  neglect  the  diffusive  contributions  to  the 
^  flux;  that  is,  the  M  term  in  the  collision  integral,  and  the  term  -(2/3)5(wg)/5w  in  the 

electric  field  expression.  The  spectrum  is  then  a  balance  between  the  mean  rates. 

The  mean  loss  rate  NqvL  increases  with  increasing  energy,  and  the  mean 
gain  rate  e*E*/mi/m  decreases  with  increasing  energy.  The  rates  cross  and  are  equal 

•  at  some  energy  Wi.  An  electron  with  energy  w  >  Wi  loses  energy  faster  than  it  gains 
it,  and  one  with  energy  w  <  wi  gains  energy  faster  than  it  loses  it.  Consequently  all 
electrons  end  up  with  energy  Wi;  the  spectrum  collapses  to  a  delta  function  at  the 
energy  where  the  rates  cross.  Thus  the  physical  process  which  actually  determines 
the  energy  spectrum  is  spread  in  gain  and  loss  rates  about  the  mean.  Mean  rates 

•  determine  only  the  average  energy  in  the  swarm,  cf  Figures  la  and  lb.  It  is  clear  that 
for  electrons  with  a  spread  in  energies  to  be  in  steady  state  there  must  be  a  balance 
between  gain  and  loss  rates  over  a  spread  in  energies,  and  this  spread  determines  the 
spectrum. 


This  shows  that  in  the  CSDA  (M  =  0),  the  only  physical  process  which 
establishes  the  spectrum,  that  is,  gives  body  to  the  distribution  function,  is  energy 
diffusion  due  to  elastic  scattering  in  the  electric  field,  since  in  the  collision  integral 
only  the  mean  loss  rate  is  retained.  The  CSDA  is  poor  because  it  neglects  “half”  of 
the  physics  determining  the  spectrum.  We  shall  later  show  it  also  violates  detailed 
balance. 


Since  energy  diffusion  is  so  important,  we  must  keep  the  diffusive  contribu¬ 
tion  in  Cg  as  well  to  deduce  a  realistic  spectrum.  While  L  is  the  lowest  order  term  in 
the  expansion  of  Cg,  M  is  the  lowest  order  contributor  to  the  energy  distribution. 
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y' 

RESULTANT  SPECTRUM 

Figure  la.  Balance  of  mean  gain  and  loss  rates.  Resultant  spectrum  is  a 
S  function  at  Wi. 


Figure  lb.  Balance  of  complete  gain  and  loss  rates.  CSDA  neglects 
spread  in  energy  loss. 
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From  this  point  of  view  X  should  not  be  considered  just  another  term  in  a 
Taylor  Series  expcinsion  of  a  function.  Rather  it  represents  a  critical  piece  of  physics 
that  must  not  be  neglected.  The  functional  form  in  which  it  appears  is  a  Fokker- 
Planck  approximation  to  the  true  random  walk.  The  approximation  is  best  when 
fractional  energy  loss  per  collision  is  small,  and  large  energy  changes  occur  only  by  a 
succession  of  small  ones. 

In  particular,  A(  allows  some  particles  to  lose  energy  slower  than  the  mean, 
thereby  supporting  the  high  energy  tail,  and  greatly  increasing  the  avalanche  rate  or 
any  process  depending  on  the  higher  energy  parts  of  the  spectrum. 

Conventional  analytic  treatments  of  the  collision  integral®’^  retain  only  the 
mean  energy  loss,  often  expressing  that  due  to  inelastic  processes  as  an  adjustable 
multiple  of  that  due  to  elastic  processes.  These  treatments  miss  an  important  process 
contributing  to  the  spectrum. 

It  is  just  as  common  for  discrete  losses  due  to,  say,  vibrational  or  electronic 
transitions,  to  not  be  included  in  a  local  differential  analysis,  it  being  believed  theit* 
discrete  nature  precludes  a  differential  treatment.  Brief  attention  has  been  paid*  to 
inelastic  collision  effects  on  the  distribution  function  when  molecular  energy  states  are 
closely  spaced,  as  for  rotations  of  diatomic  molecules,  but  a  more  general  treatment 
such  as  that  presented  here  seems  not  to  have  been  pursued.  The  all  important  role 
played  by  spread  of  energy  loss  about  the  mean  seems  not  to  have  been  appreciated. 

In  fact,  so  long  as  average  fractional  energy  loss  per  collision  is  small  com¬ 
pared  with  unity,  as  is  the  case  in  Nj  and  especially  O2  over  wide  energy  ranges,  a 
differential  treatment  including  both  L  and  M  is  justified,  and  its  accuracy  for  cal¬ 
culating  transport  coefficients  is  borne  out  by  direct  calculation.  This  would  almost 
certainly  be  true  as  well  in  many  diatomic  and  polyatomic  molecules. 


2.9  ENERGY  LOSS  PER  COLLISION 

Electrons  perform  a  random  walk  in  energy  space  as  they  gain  energy  from 
E  and  lose  it  to  inelastic  collisions.  The  mean  energy  loss  per  energy-  transfer  collision 


Wt  = 
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The  fractional  energy  loss  ^  =  w^/w  in  N*  and  Oj  at  300K  is  shown  in 
Figure  2.  ^  includes  inelastic  and  super-elastic  collisions  with  molecules  in  thermally 


10-»  10-'  1  10 
Electron  Energy  (eV) 

Figure  2.  Fractional  energy  loss  per  energy  transfer  collision  at  300K  in 
N2  and  Oj. 

populated  rotational  states.  ^  is  seen  to  be  safely  small  except  below  about  0.01  eV, 
where  rotational  transitions  with  energies  4JBo  ~  .01  eV  are  common,  and  near  10  eV 
where  high  electronic  excitations  are  most  important.  (Here  Bo  ~  2  x  10~*  eV  is  the 
rotational  constant,  and  J  10  is  a  typical  populated  rotational  quantum  number.) 
Except  for  these  energy  regions  one  would  expect  the  Fokker-Planck  approximation 
to  be  quite  good. 

Even  at  energies  large  enough  for  electronic  excitation,  which  cause  large 
single  energy  jumps,  vibrational  transitions  still  occur,  reducing  the  average  fractional 
energy  loss.  Figure  2  shows  at  10  eV,  say,  that  this  appears  to  be  more  true  for  Oj 
than  N2.  Some  computed  swarm  parameters  will  depart  from  experimental  values 
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C  in  Nj  when  energy  jumps  due  to  electronic  transitions  are  important  in  reducing  the 

tail,  or  when  the  bulk  of  the  spectrum  spans  the  2.5  eV  resonance  spike.  Thus  the 
characteristic  energy  will  be  computed  to  be  too  large  waen  it  is  near  1  or  2  eV. 
Agreement  is  better  in  Oj. 


2.10  EQUATION  FOR  TIME  DEPENDENCE  OF  ENERGY 
SPECTRUM 


Equation  20  controls  the  time  evolution  of  the  spectrum  in  time  varying 
fields  in  the  absence  of  sources  and  sinks: 


^  A 

dt  dvf  I  raUm 


2a., 

8  - 


-  NovLg  -  No 


dvf 


=  0 


(25) 


Sources  and  sinks  are  to  be  added  to  the  right  hand  side  for  a  specific  problem  as' 
appropriate.  Eq.  25  is  valid  for  E  varying  slowly  compared  with  the  momentum 
transfer  frequency  Um,  but  on  any  time  scale  relative  to  the  mean  energy  transfer 
frequency.  For  example,  following  a  pulsed  ionization  source  this  equation  determines 
the  relaxation  of  the  spectrum  toward  equilibrium  in  the  electric  field  and  warm  gas. 

Clearly  Equation  25,  with  vanishing  right-hand  side,  preserves  the  normal¬ 
ization  n  =  / gdw.  The  swarm  average  energy 


wg  dw 


(26) 


is  then  readily  shown  to  obey 


dw  ^  _  1  /  ,  .r  r 

—  =  -eE  •  Vd - /  dw  NovLg 

at  n  J 


(27) 


in  the  absence  of  sources  and  sinks.  The  swarm  theory  energy  equation  is  obtained 
by  introducing  a  swarm  energy  transfer  collision  frequency  i/w  by 


i  dw  NovLg  =  ^w  -  ^T^  i/w  , 


(28) 
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where  T  is  gas  temperature,  and  using  Eq.  (28)  to  replace  the  last  term  in  Equa¬ 
tion  (27).  The  direct  contriubtion  of  X  to  dw/dt  has  dropped  out.  At  still  affects  g 
and  therefore  the  value  of  vj  and  i/*  through  their  definitions  in  terms  of  the  spectrum, 
and  so  indirectly  alters  dvi/dt. 

For  situations  where  E  =  0,  the  equation  resulting  from  Eq.  (25), 


(29) 


determines  the  relaxation  of  g  from  an  initial  condition  toward  thermal  equilibrium. 
Without  the  E  term,  Eq.  25  is  no  longer  limited  by  Vdnft  Vthermai  (two  term  Leg¬ 
endre  series),  and  so  applies  to  electrons  of  any  energy.  It  includes  mean  energy  loss, 
straggling,  and  thermal  heating  (through  thermally  populated  rotational  states  in  L 
and  M). 


If  straggling  were  to  be  neglected  (CSDA),  the  resulting  equation 


I  -  =  0  (CSDA) 


(30) 


would  show  simple  convection  in  energy  space,  and  therefore  laminar  flow.  Electrons  of 
higher  energy  would  never  overtake  ones  of  lower  energy.  Straggling  destroys  laminar 
flow,  and  permits  passing  in  energy  space  as  actually  occurs  in  a  random  walk. 

Equation  25  can  be  written  «is  well  for  F, 


dY 

No^ 


V  d 
w  dvf 


,  w 


L  + 


1  5(wA() 


w  5w 


F  +  w 


(eE/No)^ 

3<7m 


+  At 


Sl=" 


(31) 


which  will  permit  simpler  understanding  later. 


2.11  STEADY  STATE  SPECTRUM 


When  d ld\.  =0,  the  steady  state  solution  of  (31)  is 
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F  =  Fo  exp 


-dw 


(32) 


where  Fq  is  a  normalization  constant.  This  is  an  improvement  over  the  CSDA, 

F  =  F„exp{-/;i^^dw|  (33) 

which  results  when  At  is  neglected. 

The  value  of  Eq.  (32)  lies  both  in  its  simplicity  and  (approximate)  correct¬ 
ness,  as  well  in  its  directly  exhibiting  the  dependence  of  the  distribution  function  on 
the  most  important  molecular  parameters  L,  and  At,  and  on  the  experimental  pa¬ 
rameter  E.  Scaling  laws  to  other  gases  therefore  become  more  transparent.  Likewise, 
the  dependence  of  transport  coefficients,  which  are  weighted  integrals  over  F,  on  these 
parameters  is  also  exhibited.  Section  2.12  will  further  show  the  intuitive  correctness 
ofEq.  (32). 

The  Druyvesteyn  spectrum  is  itself  a  special  case  of  Eq,  33  for  =  constant, 
and  elastic  collisions  only.  In  this  case 


L  = 


2m 


(34) 


Then  F  reduces  to 


F  =  Foexp 


3ma^w*  ] 

M(eE/No)*  J 


=  Fd 


(35) 


In  this  case,  the  more  correct  Eq.  (32)  is  essentially  the  same  as  Eq.  (35)  because,  for 
elastic  collisions,  A(  ~  i(2m/M)wL  <.  wL,  and  the  Af  correction  is  negligible.  Thus 
in  this  special  case  the  improved  spectrum  does  not  mar  already  good  agreement  with 
data. 


Later  sections  compute  swarm  parameters  with  the  spectrum  Eq.  (32)  for 
N2  and  Oj. 
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2.12  DIFFUSION-CONVECTION  THEORY 


The  functional  form  for  the  steady  state  spectrum  Eq.  32  is  easily  understood 
in  terms  of  elementary  diffusion  and  convection  theory. 

Let  a  quantity  0  diffuse  and  convect  in  one  dimension  x  with  diffusion 
coefficient  D(x)  and  convection  velocity  V(x): 


at  ax^  ax  ax 


(36) 


In  steady  state  this  equation  shows  that  the  flux  must  be  constant  in  t  and  x: 


V0  —  D-^  =  constant  =  0  ,  (37) 

ax 

where  the  constant  is  taken  to  be  0  since  we  assume  0  and  dxiijdx  vanish  at  an  end 
point.  Then 

=  X 

0ax  "  D 

In  steady  state  the  convection  velocity  and  diffusion  coefficient  determine  the  loga¬ 
rithmic  derivative  of  0  by  Eq.  38.  Consequently  0  is  distributed  according  to 

0  =  0oexp|y^Xdx|  (39) 

If  we  now  take  x  to  be  the  electron  energy  w,  and  0  the  spectrum  F,  the  flux  in 
curly  brackets  in  Equation  (31)  identifies 


V  =  -Nov 


D  =  Nov 


(eE/No)* 


3<7„ 


-l-M 


(40) 


(41) 
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so  that  the  solution  Eq.  39  is  the  steady  state  spectrum  Eq.  32. 

-  NqvL  is  the  prescribed  convection  velocity,  and  the  (1/w)  5(wA()/5w  term 
is  the  additional  convection  due  to  the  energy  dependence  of  the  diffusion  coefficient. 
Their  sum  is  the  total  convection  velocity.  The  first  term  in  Eq.  41  is  the  diffusion 
coefficient  due  to  elastic  scattering  in  E,  and  the  M  term  is  that  due  to  energy  loss 
straggling.  Their  sum  is  the  total  diffusion  coefficient. 

The  electric  field  Ed  at  which  the  two  diffusion  coefficients  are  equal  is 

eEn  /  ~ 

,  (42) 

and  is  a  function  of  w.  In  O2  it  is  as  small  as  0.1  Td  for  w  ^  0.1  eV,  about  1  Td 
near  0.3  eV,  10  Td  at  4.0  eV,  and  greater  than  100  Td  for  w^  9  eV.  In  Nj  it  is  more 
widely  variable  because  of  the  resonance  spike  which  causes  both  and  .M  to  peak 
near  2.5  eV. 

(The  mean  gain  velocity  e*E*/mi/m  and  the  velocity  -(2/3)e*E*/mi/ni  due 
to  the  energy  dependence  of  the  electric  field  diffusion  coefficient  would  also  appeau’ 
directly  if  we  had  applied  these  ideas  to  Equation  (25)  for  g.  They  drop  out  in  the 
numerator  of  the  final  distribution  function^). 


2.13  EXACT  LIMIT  FOR  SMALL  QUANTUM  TRANSITIONS 

When  E  =  0,  the  spectrum  Eq.  32  must  reduce  to  a  Maxwellian.  Reference  2 
derives  explicitly  the  condition 


L  +  _  1 

W  T 

that  assures  this.  Here  T  is  gas  temperature.  This  condition  is  exactly  true  only  for 
the  limit  of  small  molecular  quantum  transition  energies;  it  is  only  approximately  true 
for  real  cross  sections  and  real  L  and  >f.  Its  derivation  requires  detailed  balance,  and 
uses  the  complete  expression  for  L 


L(w) 


ENi(Ek-Ei)aik(w) 

i,k 


17 


(44) 


I 


I 


=  A^Die-^‘/'^(Ek-Ei)aik(w) 
>.k 


that  includes  thermally  populated  excited  states  of  energy  E;,  and  degeneracy  D,,  £ind 
with  electronic  cross  section  <7ik(w)  for  transitions  from  i  to  k.  An  expression  similar 
to  Eq.  44  holds  for  H.  Ni  =  ADie~®‘^^  is  the  population  of  the  i-th  level,  where 
A=No/E.  Die-®=>/T. 

When  |Ek  — Eil  T  or  w,  then  an  expansion  in  the  small  ratios  |Ek  —  Ei|/T 
and  |Ek  — Ei|/w,  together  with  detailed  balance,  leads  to  condition  (Eq.  43).  This  con¬ 
dition  is  the  necessary  relation  that  must  hold  between  the  Fokker-Planck  coefficients 
of  friction  and  diffusion.  In  particular,  the  CSDA  (L  7^  0,  =  0)  violates  it  and  is 

therefore  inconsistent  with  detailed  balance. 

Eq.  32  appears  to  be  significantly  superior  to  the  CSDA,  the  fractional 
energy  loss  i  being  small  enough  so  that  the  Fokker-Planck  approximation  is  valid. 
When  used  to  compute  transport  coefficients  by  integrating  over  the  spectrum,  a 
higher  accuracy  is  achieved,  and  results  are  in  surprisingly  good  agreement  with  data. 

In  this  limit,  using  Eq.  (43),  one  could  eliminate  L  in  Eq.  (32)  to  obtain 
F-Foexpl-i/',^  dwj  (45) 

I  3<7n,M  ^  ^  ) 

reducing  the  dependence  on  molecul2U’  and  experimental  quantities  to  the  single  pa¬ 
rameter  (eE/No)*/3amA(.  This  form  is  appealing  in  its  simplicity,  and  is  essentially 
the  Chapman-Cowling  result  (discussed  in  Section  2.14)  with  LT  replaced  by  AI.  How¬ 
ever,  it  focusses  attention  on  the  straggling  function  M,  £ind  detracts  from  the  more 
physically  transparent  form  of  Equations  (32)  and  (39).  We  shall  therefore  work  with 
the  distribution  function  in  the  form  of  Equation  (32). 


2.14  PREVIOUS  DISTRIBUTION  FUNCTIONS  AS  SPECIAL  CASES 


W’  have  already  mentioned  in  Section  2.11  that  the  Druyvesteyn  spectrum  is 
a  special  case  of  the  CSDA  (for  =  constant),  which  is  a  special  case  of  Equation  (32). 
Other  historically  important  distribution  functions  have  been  obtained  by  Pidduck®; 
Davydov*®;  Morse,  Allis,  and  Lamar**;  Chapman  and  Cowling**;  and  by  Wannier**. 
All  of  these  are  intended  to  apply  in  the  case  of  elastic  collisions  only.  The  later  work 
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permits  the  momentum  transfer  cross  section  crm(w)  to  be  a  function  of  energy,  and 
accounts  for  thermal  gas  motion.  We  now  show  these  are  all  approximate  special  cases 
of  Equation  (32). 

Huxley  and  Crompton  show  that  all  formulae  are  included  in  that  of  Chap¬ 
man  and  Cowling’s,  which  in  our  notation  is 


Fee  =  Foexp 


mvdv 

MfeE/No)^ 

3(inv<»„)* 


(46) 


where  T  is  gas  temperature,  and  M  is  the  molecular  mass.  Thus  we  need  only  show 
that  this  equation  is  contained  in  Equation  (32). 

For  the  conditions  stated,  L  is 


T  2m 


and,  using  mvdv  =  dw.  Equation  (46)  can  then  be  written 


(47) 


Fee  =  Foexp 


-/• 


Ldw 


lll/Noi?  +  XL 


(48) 


That  this  is  a  special  case  of  Equation  (32)  may  be  seen  by  using  Equa¬ 
tion  (43)  to  rewrite  the  At  in  the  denominator  of  the  integrand  of  Equation  (32), 


w  aw 


(49) 


so  that  it  becomes 


F  =  Foexp 


M; 


L  1  1  g(wX) 

_ w  dw 


dw 


(50) 


differing  from  Equation  (48)  only  by  the  {l/w)d{v/M)/dvf  added  to  L.  But  for  elastic 
collisions. 
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(51) 


and 


u  12m 

^  =  2  m"'' 


1 


lA 

w  dyf 


(wM) 


2m  w  d\j 

M  2L  dv! 


<C  L 


(52) 


Therefore  the  M  correction  can  be  dropped,  and  the  two  expressions  become  the  same. 

For  elastic  collisions  only,  collisional  energy  transfers  are  indeed  small  com¬ 
pared  with  w  and  T,  so  we  are  justified  in  using  Equation  (43).  Also  for  elastic 
collisions  the  M  correction  is  small,  allowing  the  above  approximations. 

When  inelastic  collisions  contribute,  as  in  vibrational  excitations  of  diatomic 
molecules,  the  contribution  {l/w)d{wM)/dw  to  the  convection  velocity  L,  arising  from 
the  energy  dependence  of  the  diffusion  coefficient,  is  not  necessarily  small,  and  can 
no  longer  be  neglected.  The  distribution  function  in  Equation  (32)  extends  earlier 
theoretical  expressions  to  the  case  of  (small)  inelastic  collisions  as  well,  and  contains 
the  earlier  work  as  special  cases. 


2.15  STEADY  STATE  RESULTS  FOR  N, 

Earlier  reports^-’  collect  requisite  cross  sections:  momentum-transfer,  ro-  ^ 

tational,  vibrational,  and  electronic.  Figure  3  shows  L  and  M  for  Nj  at  300K 
constructed  from  them,  including  thermally  populated  rotational  states  and  super¬ 
elastic  rotational  collisions.  The  resonance  spike  near  2.5  eV  and  the  high  first  elec¬ 
tronic  excitation  (6.17  eV)  make  Nj  a  particularly  stressful  case  in  which  to  test  the 
Fokker-Planck  Equation  32.  In  fact  the  spectrum  computed  in  a  field  of  100  Td,  Fig-  ® 

ure  4  [Reference  1), while  superior  to  the  CSDA,  is  not  very  close  to  more  accurate 
published  calculations**.  (ITd  =  10“*^V  —  cm*).  However  transport  coefficients  are 
in  much  better  agreement  with  data,  except  for  those  dependent  on  the  tail. 

Figure  5  shows  the  calculated  drift  velocity  (open  circles)  with  data  compiled  • 

by  Dutton*®  over  some  five  decades  in  E/Nq.  The  agreement  is  excellent  everywhere 
except  between  10"*®  and  10"*®Vcm*,  where  calculated  drift  velocities  are  some  20% 
low.  At  these  field  strengths,  the  characteristic  energy  is  between  1  and  3  eV,  and 
is  computed  to  be  too  large  (Figure  6)  partly  because  of  the  high  spectral  tail,  and 
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w  (eV) 

Figure  3.  At  and  L  functions  for  N2  at  300  K. 

partly  because  the  spectrum  straddles  the  resonance.  At  even  higher  fields  the  drift 
velocity  returns  to  good  agreement. 

Figure  7  shows  swarm  momentum-  and  energy-transfer  collision  frequencies 
as  a  function  of  characteristic  energy  up  to  1  eV,  compared  with  more  detailed  calcu¬ 
lations.  Again  the  agreement  is  excellent  except  for  the  energy-transfer  frequency  at 
1  eV. 

Figure  8  shows  the  avalanche  rate  above  10“^®V  —  cm*.  At  relatively  low 
E/No  it  is  computed  to  be  more  than  an  order  of  magnitude  too  large,  due  again  to 
the  high  tail  seen  in  Figure  4. 
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F  (eV-^/2) 


Figure  4. 


w{eV) 


Energy  spectrum  at  10  V-cm*  from  Eq.  (32),  the  CSDA, 
and  Reference  14. 
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E/N  (10*'’V  ai*) 

Figure  5.  Calculated  Nj  drift  velocity  (open  circles)  with  data  compiled 
by  Dutton  (Kef.  15). 


0,/«  (w«|ts)  D/tf  (w«lu) 


(V/CM  t*fp) 

10**  I  O'*  to*' 


Figure  6.  Calculated  gwarm  characteristic  energies  in  Nj  (open  circles) 
with  data  compiled  by  Dutton  (Ref.  15).  Ambient  temperature 
(300  K). 
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Figure  7.  Swarm  momentum  transfer  (u„/N)  and  energy  transfer  (t^u/N) 
collision  frequencies.  The  large  open  circles  are  our  calculated 
values,  superimposed  on  a  graph  from  Frost  and  Phelps  (Ref. 
16).  The  solid  curves  are  averaged  experimental  data  at  300  K. 
The  other  four  sets  of  points  are  various  theories  within  Ref. 
16. 
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2.16  STEADY  STATE  RESULTS  FOR  Oj 

O2  cross  sections  were  collected  in  previous  reports^’^.  There  is  a  discrepancy 
in  the  momentum-transfer  cross  section,  it  being  more  poorly  understood  in  Oj  than 
N2.  We  choose  two  alternate  sets  as  in  Reference  2.  The  first  is  the  set  recently 
recommended  by  Phelps^^,  denoted  in  Figures  9  through  14  by  open  triangles;  and 
the  second  is  that  of  VVadzinski  and  Jasperse^*  corrected  by  Kieffer’s^®  low  energy  tail, 
denoted  by  open  circles. 

The  rotational  cross  sections  are  those  of  Gerjuoy-Stein  with  quadrupole 
moment  Q  =  1.8.  For  the  vibrational  cross  sections,  sharp  spikes  are  smoothed  out, 
choosing  their  energy  integal  to  agree  with  Phelps’  normalization.^^  The  electronic 
cross  sections  are  those  compiled  earlier  by  Archer*®  zmd  in  use  by  the  DNA  auroral 
physics  community  for  years. 

The  resulting  L  and  At  are  shown  in  Figure  9  at  300K.  Rotational  transitions 
are  responsible  for  almost  all  energy  loss  below  about  0.2  eV.  The  hump  between  0.2 
and  1  eV  is  due  to  vibrational  transitions. 

The  spectrum  F  at  E/Nq  =  10  and  100  Td  is  shown  in  Figure  10,  along 
with  that  from  the  CSDA.  Serious  differences  between  the  two  are  apparent  at  higher 
energies.  Since  many  transport  coefficients  do  not  depend  sensitively  on  the  high 
energy  tail  the  CSDA  is  not  wildly  different  from  either  the  L,  A(  approximation  or 
data,  except  for  avalanching  at  low  values  of  E/Nq. 

Attachment  has,  of  course,  been  ignored  for  the  purpose  of  computing  F, 
and  we  have  found  no  published  O2  spectra  computed  with  the  full  collision  integral 
with  which  to  compare  the  L,  .M  approximation. 

Figures  11,  and  12  show  the  characteristic  energy  zmd  drift  velocity,  drawn 
as  circles  and  triangles  over  data  compiled  by  Dutton.  Below  about  10“^*Vcm*,  dif¬ 
ferences  in  drift  velocity  between  the  two  sets  of  momentum  transfer  cross  sections, 
with  each  set  agreeing  with  one  set  of  data,  reflect  difficulties  in  performing  drift  ve¬ 
locity  measurements  in  the  electro-negative  gas.  Agreement  is  excellent  over  some  five 
orders  of  magnitude  in  E/Nq. 

Figure  13  is  the  momentum-  and  energy-transfer  collision  frequencies  shown 
with  more  accurate  calculations  by  Hake  and  Phelps**.  With  the  momentum  transfer 
cross  section  of  Wadzinski  and  Jasperse**’**  the  agreement  is  striking,  except  in  the 
energy-transfer  frequency  at  5  eV  where  electronic  transitions  cause  a  larger  energy 
transfer  rate  than  the  Fokker-Planck  expression  predicts. 


Electron  energy  (eV) 

Figure  10.  Energy  spectrum  in  Oj  at  E/Nq  =  10  and  100  Td. 
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SECTION  3 


TIME  DEPENDENT  SPECTRA 


Equation  (31),  together  with  sources  or  sinks  on  the  right  hand  side,  de¬ 
scribes  the  time  evolution  of  the  distribution  function  F(w,t)  in  a  constant  or  time 
varying  electric  field.  It  must  be  solved  numerically.  In  this  section  we  record  the 
numerical  integration  scheme  we  have  found  most  successful. 

We  need  to  solve 


^  -  A 

dr  ^  dw 


qF  -h  r 


A' 

dw 


(53) 


where 


T  =  Not 


(54) 


p  =  v/w 


(55) 


q  =  w 


L  +  ~~(viM) 
w  aw 


(56) 


r  =  w 


(eE/No)» 


+  M 


(57) 


with  the  initial  condition  F(w,0)  =  Fo(w),  over  an  energy  range  from  w=0  to  w=Wniax- 
Wmax  depends  on  the  problem,  especially  the  sources,  but  it  will  be  a  few  eV  or  more. 


3.1  STABILITY  CONDITIONS 


Equation  (53)  has  an  advective  term  and  a  diffusion  term.  For  explicit 
schemes  a  Courant  stability  condition 


Aw 
At  < - 

pq 


(58) 


and  a  diffusion  stability  condition 


At  <  - —  (59) 

2pr  ^  ^ 

have  to  be  met.  An  unconditionally  stable  implicit  difference  scheme  would  relax  one 
or  both  stability  criteria,  but  the  accuracy  condition  is  about  the  same;  it  would  be 
dangerous  to  exceed  these  criteria  by  too  much. 

From  0  to  Wmax  Figures  3  and  9  show  q  and  r  vary  over  several  orders  of 
magnitude.  If  one  were  to  difference  in  w  with  constant  energy  spacing  Aw,  there 
would  result  an  unacceptably  small  Courant  time  step  at  larger  w  with  the  large 
velocity  pq  and  diffusion  coefficient  pr.  Choosing  increasing  Aw  would  help  solve 
this,  but  at  the  expense  of  lost  centering  in  the  spatial  derivatives  and  therefore  only 
first  order  accuracy. 


3.2  CHOICE  OF  INDEPENDENT  VARIABLE 


Both  problems  can  be  alleviated  by  using  the  independent  variable 


u  =  £n(w  +  Wo)  , 

where  wq  is  a  constant,  so  that  Equation  (53)  becomes 


aF  d  f  d¥\ 


where 


(60) 


(61) 


P  =  pu' 


Q  =  q 


(62) 

(63) 
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R  =  ru' 


(64) 


u'  =  du/dw  = 


w  +  Wo 


=  e 


(65) 


The  stability  criteria  become 


and 


At  < 


Au 

PQ 


Au 

pqu' 


(w  +  wq)  Au 

pq 


At  < 


Au* 

2PR 


Au*  _  (w  +  wo)*Au* 
pru'*  pr 


(66) 


(67) 


effectively  making  Aw  proportional  to  w  when  one  chooses  constant  Au.  wo  is  chosen 
larger  than  0  so  that  w  =  0  is  a  computation  point.  Differencing  Equation  (61)  with 
constant  Au  allows  both  spatial  centering  for  second  order  accuracy  and  effectively 
increased  cell  size  for  larger  time  steps. 


3.3  PHYSICAL  TIME  SCALES 


The  time  scale  for  an  electron  to  change  its  energy  by  5u  is  of  order 


T  ~  5u/PQ  ~  5w/pq  (68) 

in  the  first  term,  and 

T  ~  6u*/PR  ~  5w*/pr  (69) 

in  the  second  term.  These  time  scales  must  be  resolved  for  an  accurate  solution. 
But  both  of  them  are  long  compared  to  the  time  steps  in  Equations  (66)  and  (67) 
if  Au  <  6u.  Consequently,  simply  choosing  Au  to  resolve  energy  scales  of  interest 
guarantees  a  At  that  assures  adequate  resolution  in  the  time  evolution  of  F  as  well. 
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3.4  NORMALIZATION 


Equation  (53)  preserves 


fOO 

O  ^  ~  J  V^F(w)dw  (70) 

if  there  are  no  sources  or  sinks  on  the  right  hand  side,  and  it  is  important  that  the 
difference  algorithm  does  also.  Experience  with  non-conservative  difference  schemes 
has  shown  that  they  do  not  preserve  rj  accurately  enough.  We  shall  therefore  choose 
a  flux  conservative  differencing  scheme. 


3.5  DIFFERENCE  SCHEME 


Define  energies  Wj  at  points  i  =  1,2,- 


•,N  by 


Wi 

=  0 

• 

Ui 

=  £n(wi  -1-  Wo)  =  £n(wo) 

Ui 

=  ui -1- (i  -  l)Au  ,  i  =  2,3,  •••,N 

Wi 

II 

C 

1 

0 

(71) 


•  with  un  =  fn(wmax  +  Wo),  Au  =  (un  —  ui)/(N  —  1).  Define  F  at  the  same  points, 

Fi  =  F(uO  . 

•  Integrating  P“^  times  Equation  (61)  between  Uj- 1  Au  =  Ui-i/j  to  Ui  -|- 1  Au  =  Ui+1/2 
gives  the  basic  equation 


QF-i-R 


d\i 


i+1/2 

i-1/2 


(72) 


where 
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r'+i/s  du  /*+!/*  w  , 

Ayi  =  /  -5-  =  /  -dw  , 

Ji-l/S  P  Ji-l/l  V 


and  we  approximate  the  normalization  element  by 


Ji-l/2  P 

The  form  of  Eq.  (72)  is  flux  conservative,  preserving 


(73) 


(74) 


r  rw 

I  v^wEdwoc  /  — Fdw  =  ^  AyjFi  +  boundary  terms  , 

^  i=2 

except  for  fluxes  controlled  by  boundary  conditions,  since  Eq.  (72)  says 


(75) 


a 

Ir  ?  = 


df 


N-1/2 


l+i/J 


(76) 


ting 


The  Crank-Nicholson  implicit  difference  algorithm  is  used  for  Eq.  (72).  Set- 


F"  =  F(wi,r„)  , 


(77) 


Eq.  (72)  is  advanced  from  r„  to  r^+i  =  Tn  +  Ar  by 

_  Fl'^^  -  FP 
dr  At 


(78) 


and  the  two  terms  on  the  right  hand  side  at  i  +  1/2  are  averaged  at  n  and  n+1. 


(QF) 


fFr+Fr„  ,  Fr‘+FKi‘) 

i+i/i  -  jQi+i/i  I  — 2 —  — 2 — )  ’ 


(79) 


38 


R 


i+l/2 


i+l/2 


Fr+1 


Au 


T?n  ipn+l  T7 

Il  +  rh+x  Fiil  -F 


■i+l/2 


Au 


n+l  'X 


The  right  hauid  side  at  i  -  1/2  is  the  same  with  i  replaced  by  i  -  1. 
Boundary  conditions  are  that  the  flux  vanishes, 


(80) 


dF 

QF  +  R— 
ou 


1+1/2 


—  Qi+1/2 


pn+l  ^  pn+1 


I  pn+1 
^1  +  1/2 


pn+1  _  pn+1 


Au 


=  0 


(81) 


dF 


N-1/2 


Fn+1  ■  pn+1  pn+1  pn+1 

N-l  +  ^N  ,  r>n+l  ~^N-1 

=  V«N-l/2  ^  +  Rn_i/, - - 


=  0  ,  (82) 


which  guarantees  normalization  preservation. 
The  resulting  implicit  equations 


aiF'VY  +  AFf^'  +  7iFr_V  =  Zi  ,  i  =  1, 2,  •  •  • ,  N 


(83) 


r^^ire  inverting  only  a  tridiagonal  matrix.  The  coefficients  a\,P\,  and  7,  are  readily 
solved  for  in  terms  of  the  Q’s  and  R’s,  and  Z\  involves  only  quantities  at  time 
Boundary  conditions  Eq.  (81)  and  (82)  are  incorporated  in  Eq.  (83).  7i  =  on  =  0. 

The  first  nonzero  computation  energy  is 


wj  =  e“*  —  Wo  =  —  Wo  =  Wo(e^“  -  l) 

«  WoAu  (84) 

since  Au  «C  1.  With  N~  10*,  and  w^ax  a  few  eV,  Au  is  typically  Au  ~  .05.  Then  Wo 
essentially  controls  Wj.  Desiring  Wj  ~  10“*  eV  means  wo  should  be  chosen  of  order 
.02. 


For  the  first  few  cells  iAu  •<  1,  and  Wj  is  equally  spaced  with  spacing 
Aw  ~  Wj.  Then  the  energy  spacing  starts  to  increase,  with 

I 
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Awi  =Wi+i-Wi  =  =e“‘(e^“-l)  =  (wi  +  Wo)(e^“ - 1)  «  (wi+wo)Au  (85) 

approximately  linearly  increasing  with  w  when  Wj  >  wq. 

We  now  investigate  time  dependent  spectra  with  and  without  an  electric 
field  using  numerical  solutions  of  Eq.  (31). 


SECTION  4 


RELAXATION  OF  THE  SPECTRUM  IN  THE  ABSENCE  OF  FIELDS 


In  the  absence  of  an  electric  field,  collisions  cause  an  initial  distribution  to 
relax  toward  thermal  equilibrium  in  the  host  gas,  a  process  described  in  the  L,  At 
approximation  by  Equation  (29),  or  Equation  (31)  with  E  =  0.  The  spectrum  starts 
at  F(w,t=0)  =  Fo(w).  We  can  study  in  particular  the  difference  in  the  evolving  spectra 
between  the  CSDA  and  when  the  correct  At  is  used. 

We  use  the  numerical  solution  of  Equation  (53)  described  in  Section  3.  With 
E  =  0  the  formulation  is  not  restricted  to  the  first  two  Legendre  terms,  and  so  applies 
to  electrons  of  any  energy,  so  long  as  their  fractional  energy  loss  per  collision  remains 
small.  Nevertheless  we  are  mostly  concerned  with  electrons  of  relatively  small  energies, 
the  so-called  swarm  regime.  Only  the  motion  of  primary  electrons  is  treated;  increasing 
numbers  by  avalanching  is  not  included. 

4.1  RELAXATION  IN  Oj 

Consider  a  hot  electron  swarm  in  room  temperature  O2.  We  arbitrarily 
choose  Fo(w)  to  be  a  Maxwellian  at  T,  =  0.2  eV,  and  follow  the  evolution  of  the 
spectrum  toward  a  Maxwellian  at  300K  =  0.026  eV,  ignoring  attachment. 

Figure  15  shows  the  spectrum  at  six  times,  starting  at  t  =  0.  It  begins 
to  pile  up  near  0.2  eV  because  only  relatively  slow  rotational  transitions  can  reduce 
it  from  there.  The  bunching  follows  the  small  values  of  L  seen  in  Figure  9.  It  has 
essentially  reached  its  final  Maxwellian  value  by  Not  =  2  x  10*'s/cm^. 

At  the  final  time,  between  about  0.01  eV  and  0.25  eV,  the  calculated  distri¬ 
bution  is  exp(-w/Te),  with  T,  =  .026  eV,  as  it  must  be.  At  smaller  energies  it  departs 
from  Maxwellian  because  rotational  transitions  are  comparable  to  w,  as  discussed  in 
Section  2.9.  At  higher  energies  it  is  smaller  than  Maxwellian  because  we  have  not 
included  thermally  populated  vibrational  levels  in  the  gas. 

Figure  16  shows  the  same  times  for  the  CSDA.  Comparing  with  Figure  15 
the  effect  of  M,  in  elevating  the  tail,  and  in  spreading  out  the  peaks,  is  obvious. 

A  blow  up  of  the  low  energy  region  is  shown  in  Figure  17  for  both  the  CSDA 
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Relaxation  of  the  energy  spectrum  in  O2  in  the  L,  M 
approximation.  4.4(10)  means  4.4xl0‘°s/cin*,  etc.. 
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Figure  16.  Relaxation  of  the  energy  spectrum  in  O2  in  the  CSDA.  Same 
successive  times  as  Figure  15. 
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(dashed  lines)  and  when  the  correct  M  is  employed. 

The  CSDA  spectrum  has  not  quite  collapsed  to  a  delta  function  at  the  last 
time  because  a  small  diffusion  coefficient  was  retained  for  code  stability. 

The  decay  of  the  instantaneous  characteristic  energy  e  =  Dx/m  where  Dx  is 
the  transverse  diffusion  coefficient  and  n  the  mobility,  and  the  average  energy  w,  are 
shown  in  Figures  18a  and  b.  w  differs  only  slightly  whether  is  retained  or  not,  in 
accordance  with  the  discussion  following  Equation  (34).  c  is  seen  to  differ  considerably 
more. 


The  growth  of  the  swarm  mobility  in  the  two  cases  is  shown  in  Figure  19.  A 
gas  density  of  No  =  3.36  x  10'®cm~*  heis  been  assumed  to  obtain  mobility  in  m/s/V/m. 
In  the  CSDA  the  mobility  always  lags  behind  the  correct  value.  This  result,  together 
with  a  similar  one  for  air  discussed  below,  suggests  that  mobility  or  conductivity 
calculations  in  actual  cases  may  likewise  be  inaccurate  by  comparable  amounts  if  the 
CSDA  were  used. 


4.2  RELAXATION  IN  AIR 


Figure  20  shows  the  Loss  Function  L  for  air  (.79N2  +  .21 02).  Spectrum 
rela\ation  in  the  L,  X  approximation  is  shown  in  Figure  21  with  the  complete  X, 
and  in  Figure  22  in  the  CSDA.  Differences  between  them  are  similar  to  the  case 
of  O2.  In  the  CSDA  the  distribution  above  1.5  eV  is  completely  gone  after  only 
Not  =  8  X  10®sec/cm®,  and  is  gone  above  0.5  eV  after  only  4  x  10^°sec/cm®.  The  broad 
hump  in  the  distributions  in  Figure  21  near  4  eV  is  electrons  lingering  in  the  hollow 
in  L  between  the  resonance  peak  and  the  onset  of  electronic  transitions. 


Figures  23a  and  b  show  the  decay  of  the  charsu:teristic  energy  and  the  av¬ 
erage  energy,  while  Figure  24  depicts  the  growth  of  the  mobility.  Even  though  the  L, 
X  approximation  results  in  a  distribution  with  a  higher  tail,  the  characteristic  energy 
is  dominated  by  low  energjy  electrons  and  is  higher  in  the  CSDA  (by  almost  a  factor 
of  two  at  2  xl0^*s/cm’),  causing  mobility  to  be  lower  in  the  CSDA. 
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Figure  18.  Evolution  of  the  characteristic  energy  (a)  and  average  energy 
(b)  in  Oj  in  the  L,  M  approximation  (solid)  and  CSDA 
(dashed). 
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Evolution  of  the  mobility  in  air.  L,  M  approximation  (solid); 
CSDA  (dashed). 


SECTION  5 


SPECTRAL  BEHAVIOR  FOLLOWING  A  TIME  VARYING 

ELECTRIC  FIELD 


Except  for  spectral  lag  effects  which  require  a  time  of  order  l/i^uj  the  time 
scale  for  the  drift  velocity  of  a  swarm  to  readjust  in  a  varying  electric  field  is  Ijum, 
arising  from  the  randomizing  of  directed  velocities  due  to  momentum  transfer  scat¬ 
terings.  The  time  scale  for  the  energy  spectrum  to  equilibrate  in  a  changing  field  is 
1/i^u,  the  inverse  of  the  swarm  energy  transfer  collision  frequency.  These  frequencies 
were  shown  as  a  function  of  swarm  chjiracteristic  energy  in  Figures  7  and  13  for  N2 
and  O2,  respectively. 

Generally  !/„,>>  i^u-  Thus  there  is  an  interesting  time  region  in  which 
electric  fields  vary  slowly  compared  with  Vm  but  not  slowly  compared  with  i/y.  The 
L,  A1  approximation  is  valid  in  just  this  regime,  and  can  exhibit  time  lag  effects  in 
relatively  rapidly  varying  fields. 

As  an  example,  the  EMP  from  high  altitude  nuclear  detonations  is  produced 
by  7  rays  depositing  at  altitudes  between  about  20  and  40  km  creating  Compton 
electrons  which  then  produce  secondairy  ionization.  The  secondary  electrons  are  a 
swarm,  which  drifts  in  the  self-consistent  electric  field  making  a  conduction  current. 
At  30  km,  for  example,  No  =  3.8xl0^^cm“®,  and  !/«/„  «  1.3x10“*  sec  for  characteristic 
energy  e  =  1  eV.  This  time  is  comparable  to  the  time  scale  of  the  changing  electric 
field,  and  so  the  swarm  energy  spectrum  lags  the  field  in  equilibrating  to  its  new  value. 
Thus  the  swarm  mobility  would  not  be  given  correctly  by  its  steady  state  value  in  the 
instantaneous  field. 

To  demonstrate  spectral  time  lag  effects  we  calculate  here  the  time  depen¬ 
dent  spectrum  and  transport  coefficients  of  an  originally  thermal  swarm  in  a  rapidly 
varying  pulsed  electric  field.  We  do  so  in  air  in  spite  of  the  fact  that  it  stresses 
the  Fokker-Planck  approximation,  as  mentioned  in  the  discussion  of  N2.  The  field  is 
turned  on  at  t  =  0  and  then  taken  to  vary  as  a  sine  wave  for  half  a  period: 

E  =  Eosin(7rt/To),  t  <  To 

=  0,  t  >  To  (86) 
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An  air  density  of  3.83  xlO'^cm"*  is  chosen  to  correspond  to  an  altitude  of 
30  km.  We  choose  a  peak  field  of  Eq  =  3.83  x  10®  V/m,  corresponding  to  Eo/N  =  10 
Td.  In  the  peak  field  in  steady  state  this  would  give  the  swarm  a  characteristic  energy 
of  c  =  1  eV.  Thus  a  useful  time  scale  is  l/t^u  (1  eV).  From  Figures  7  and  13,  at  the 
chosen  density,  this  time  is  «  12  ns. 

Since  the  swarm  starts  at  thermal  energies,  the  calculation  is  not  meant  to 
apply  to  high  altitude  nuclear  EMP,  but  is  intended  to  illustrate  time  lag  effects,  and 
the  usefulness  of  the  L,  At  approximation  in  quickly  calculating  spectral  behavior. 

5.1  RAPIDLY  VARYING  FIELD 

We  first  choose  To  =  2.5  x  10“®  sec,  so  the  field  peaks  in  about  l/^'u  (1  eV). 
However,  this  time  variation  is  rapid  enough  that  the  characteristic  energy,  which 
starts  at  0.026  eV,  never  rises  to  1  eV  before  the  pulse  is  over,  since  is  much  less 
for  lower  energies. 

The  field,  characteristic  energy,  and  mobility  are  shown  in  Figure  25.  e  rises 
to  only  about  0.13  eV  at  20  ns.  The  lower  dot-dashed  curve  shows  what  the  energy  ejngt 
would  be  if  the  spectrum  equilibrated  immediately  to  the  instantaneous  field.  When 
the  pulse  terminates,  e  «  0.12  eV.  At  this  energy,  l/i/^  as  150  ns,  and  the  spectrum 
and  e  subsequently  relajc  very  slowly,  lingering  for  a  long  time  at  almost  four  times  its 
final  value. 

The  mobility  /x  in  very  weak  fields  is  about  100  (m®/Vs).  It  likewise  never 
reaches  its  steady  state  value  in  the  peak  field,  achieving  a  minimum  of  only  about 
25  (m*/Vs).  The  upper  dot-dashed  curve  shows  what  the  mobility  ^ingt  would  be  if 
the  spectrum  equilibrated  instantaneously.  For  both  the  mobility  and  characteristic 
energy  there  is  almost  an  order  of  magnitude  discrepancy  during  most  of  the  pulse. 

5.2  MORE  SLOWLY  VARYING  FIELD 

As  another  example,  we  slow  the  field  down  by  a  factor  of  4,  choosing 
To  =  1  X  10"'^  sec.  This  will  allow  transport  coefficients  to  more  nearly  approach  their 
steady  state  values. 

Figure  26  shows  the  same  quantities.  The  characteristic  energy  in  this  case 
reaches  0.65  eV  at  75  ns  and  then  decays  at  a  varying  rate;  1/^  decreases  eis  e  decreases. 


Characteristic 
in  air  lasting  25  us 


The  mobility  almost  achieves  its  steady  state  value  (6  m*/Vs)  ii.  the  peak  field  before 
the  pulse  is  over.  Both  e  and  n  exhibit  a  time  lag  of  about  50  ns  relative  to  the  field, 
being  an  average  of  l/i/u  in  the  swarm  as  it  heats. 

Figure  27  shows  the  energy  spectrum  for  this  case  at  20,  40,  60,  100,  and  200 
ns,  stcirting  at  a  300K  Maxwellian  at  t  =  0  (dashed  line).  At  each  time  the  spectrum 
is  far  from  its  equilibrium  value  in  the  instantaneous  field  (not  shown). 

5.3  FINAL  COMMENT 

The  L,  A(  approximation  reduces  the  integro-differential  Boltzmann  Equa¬ 
tion  to  a  differential  equation,  which  may  be  solved  by  standard  finite  difference 
numerical  methods.  The  equation  contains  most  of  the  physics  determining  the  spec¬ 
trum,  and  is  accurate  when  the  fractional  energy  loss  per  energy-transfer  collision  is 
small.  It  is  more  accurate  than  the  CSDA,  is  conceptually  more  complete,  and  does 
not  violate  detailed  balance. 

Here  we  have  explored  only  the  cases  of  spectrum  relaxation  in  the  absence 
of  fields,  and  spectral  lag  in  rapidly  varying  fields.  One  could  easily  include  prescribed 
sources  and  chemistry  on  the  right  hand  side  to  treat  more  general  problems.  Refer¬ 
ences  1  and  2  presented  a  more  detailed  treatment  of  the  theory,  and  a  more  thorough 
discussion  of  the  steady  state  spectra  and  transport  coefficients  in  Nj  and  Oj. 
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